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ABSTRACT 

We present simulations of the compact galaxy group Stephan's Quintet (SQ) including 
magnetic fields, performed with the N-body/smoothed particle hydrodynamics (SPH) 
code Gadget. The simulations include radiative cooling, star formation and super- 
nova feedback. Magnet ohydro dynamics (MHD) is implemented using the standard 
smoothed particle magnet ohydro dynamics (SPMHD) method. We adapt two different 
initial models for SQ based on Renaud et al. and Hwang et al., both including four 
galaxies (NGC 7319, NGC 7320c, NGC 7318a and NGC 7318b). Additionally, the 
galaxies are embedded in a magnetized, low density intergalactic medium (IGM). The 
ambient IGM has an initial magnetic field of 10~^ G and the four progenitor discs have 
initial magnetic fields of 10~^ — 10"'' G. We investigate the morphology, regions of star 
formation, temperature. X-ray emission, magnetic field structure and radio emission 
within the two different SQ models. In general, the enhancement and propagation of 
the studied gaseous properties (temperature. X-ray emission, magnetic field strength 
and synchrotron intensity) is more efficient for the SQ model based on Renaud et al., 
whose galaxies are more massive, whereas the less massive SQ model based on Hwang 
et al. shows generally similar effects but with smaller efficiency. We show that the 
large shock found in observations of SQ is most likely the result of a collision of the 
galaxy NGC 7318b with the IGM. This large group-wide shock is clearly visible in the 
X-ray emission and synchrotron intensity within the simulations of both SQ models. 
The order of magnitude of the observed synchrotron emission within the shock front 
is slightly better reproduced by the SQ model based on Renaud et al., whereas the 
distribution and structure of the synchrotron emission is better reproduced by the SQ 
model based on Hwang et al.. 

Key words: methods: numerical - galaxies: spiral - galaxies: magnetic fields - galaxies: 
interactions - galaxies: kinematics and dynamics 



Stephan's Quintet (Stephan 1877; hereafter SQ), also known 
as Hickson Compact Group 92 (Hickson 1982), is the first 
discovered compact galaxy group. It is located in the con- 
stellation Pegasus. SQ consists of five galaxies (NGC 7319, 
NGC 7318a, NGC 7318b, NGC 7317, and NGC 7320, cf. 
also Fig. 1) with an estimated distance to earth of ^ 94 
Mpc (Moles et al. 1998; Appleton et al. 2006). SQ is famous 
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1 INTRODUCTION 



for a physical adjacency between four of the galaxies (NGC 
7319, 7318a, 7318b, and 7317), whereby strong interactions 
between three of these members are apparently causing tidal 
tails, a strong group- wide shock visible in X-ray (Pietsch et 
al. 1997; Sulentic et al. 2001) and radio emission (Allen & 
Hartsuiker 1972; Xu et al. 2003), and a region of active star 
formation northern of the colliding galaxies (Xu et al. 2005). 
The fifth galaxy, NGC 7320, is observed to be a much closer 
foreground galaxy and is not part of the interacting group 
(Shostak et al. 1974; Allen & Hartsuiker 1980; Moles et al. 
1997). It has a recessional velocity of ^ 740 km s~^ (Falco 
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et al. 1999). Three of the four physically related galaxies, 
the main galaxy NGC 7319 as well as NGC 7317 and NGC 
7318a, have a similar recessional velocity of ~ 6640-6670 
km s~^ (Fedotov et al. 2011) and represent the core of the 
compact group (c.f. Fig. 1). NGC 7318b is observed to be 
a high-speed intruder 5770 km s~^) and seems to inter- 
act with the IGM within the main system for the first time 
(Moles et al. 1997; Fedotov et al. 2011). Finally, there is a 
sixth galaxy, NGC 7320c, which shows a similar recessional 
velocity (~ 5990 km s~^) as the core of the compact group 
(Fedotov et al. 2011) and is therefore suggested to interact 
with NGC 7319 (Moles et al. 1997), because it reveals con- 
nected tidal features in the eastern large tidal tail (cf. the 
outer tail in Fig. 1). Therefore, it is most likely also part of 
the compact group (Arp 1973). 

Radio observations reveal the presence of magnetic fields 
in most late- type galaxies of the local Universe. The field 
strengths range from a few /xG in isolated quiet galaxies 
up to about 50 /xG in starburst galaxies (see e.g. Hummel 
1986; Klein et al. 1988; Fitt k Alexander 1993; Beck et al. 
1996; Chyzy et al. 2003). The intergalactic magnetic field 
is usually estimated to be less than 0.01 /iC (e.g. Kronberg 
et al. 2008). Galaxy interactions can cause magnetic fields 
which are much stronger compared to individual galaxies 
(Beck 2005; Drzazga et al. 2011). In this context, the SQ is 
particularly interesting because it shows both galaxy-galaxy 
and galaxy-IGM interactions. The latter is mainly visible on 
the basis of the prominent ridge of X-ray and radio emission 
crossing the system. Xu et al. (2003) estimate the minimum- 
energy magnetic field strength within the shock region of SQ 
to ^ 10 fiG. 

The morphology of an interacting galaxy system strongly 
depends on the initial properties of the progenitors. Numer- 
ical simulations can provide insights into the properties of 
the initial progenitor galaxies, as well as possible formation 
scenarios. So far, simulations of galaxy interactions were pre- 
dominantly focusing on stellar dynamics, gas fiows, star for- 
mation, supermassive black holes and feedback from stars 
and black holes (e.g. Di Matteo et al. 2005; Springel et al. 
2005a,b; Robertson et al. 2006; Cox et al. 2008; Johansson 
et al. 2009). However, over the past few years, the interest in 
simulations of magnetic fields in galactic environments has 
grown substantially (e.g. Kotarba et al. 2009, 2011; Dubois 
& Teyssier 2010; Geng et al. 2012). Synthetic radio maps 
of the interacting system of the Antennae galaxies were re- 
cently presented by Kotarba et al. (2010). They showed ra- 
dio and polarization maps calculated from their SPMHD 
simulations to be in good morphological and quantitative 
agreement with observations, implying that the simulations 
have the capability to follow the magnetic field evolution in 
a highly interacting nonlinear environment. 
Despite an enormous number of observational studies of SQ, 
revealing different features at different wavelengths, numer- 
ical simulations of this system are very rare, due to the diffi- 
culties of many free parameters in modeling such a complex 
interacting system. The simulated models for SQ commonly 
exclude the foreground galaxy NGC 7320 but include NGC 
7320c as a member of the galaxy group. Furthermore, NGC 
7317 does not show any visible features of an ongoing or past 
interaction. Therefore, it is usually not taken into account 
in the simulations. As a first basis for more detailed studies, 
Renaud et al. (2010) performed collisionless gravitational 
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Figure 1. Schematic view of the main morphological features 
of SQ. Historically, SQ consists of the galaxies NGC 7319, NGC 
7318a/b, NGC 7320 and NGC 7317. Today, however, only the 
galaxies NGC 7319, NGC 7318a/b, NGC 7317 and additionally 
NGC 7320c are observed to form an interdependent galaxy group. 
As indicated in the upper right (grey), north is up and east is to 
the left. 

N-body simulations proposing a possible formation scenario 
for SQ. Their simulation is mainly focusing on reproducing 
the stellar large-scale structure of SQ. Recently, Hwang et 
al. (2012) presented an extended model including a gaseous 
component and different galaxy models, suggesting a differ- 
ent formation history. The global morphology of the system 
is well represented in their simulation, supporting also the 
idea that the large-scale shock within the system is caused 
by the interaction of NGC 7318b with the intragroup IGM. 
However, until now, detailed numerical studies of the SQ 
concerning star formation, temperature, magnetic fields, X- 
ray and radio emission are still missing. Therefore, within 
the presented work, we perform SPMHD studies of SQ on 
the basis of the two existing models by Renaud et al. (2010) 
and Hwang et al. (2012). Thereby, we place a particular 
focus on the properties of the gaseous component, i.e. on 
the star formation rate, the temperature, and the magnetic 
field. We also investigate the presence of the shock front in 
the X-ray and radio emission. Our synthetic radio maps are 
calculated for different frequencies, allowing a better com- 
parability with observations. 

The paper is organized as follows. In section 2 we introduce 
the main observational properties of SQ. Section 3 briefiy 
summerizes the existing models of SQ and presents our ex- 
tended initial conditions. The numerical method is described 
in section 4. In section 5 we present the results of our simu- 
lated SQ models, including star formation regions, temper- 
ature distribution. X-ray emission, magnetic field structure, 
radio emission and polarization maps. We discuss the two 
different SQ models in section 6. Finally, we summerize and 
discuss our findings in section 7. 



2 OBSERVATIONAL PROPERTIES OF THE 
INTERACTING GROUP SQ 

SQ is probably the most studied compact galaxy group. The 
system shows interesting tidal features resulting from one or 
more strong interactions of the galaxies, namely the outer 
tail and the inner tail (cf. Fig. 1). It was suggested that the 
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Figure 2. R-band CCD image overlaid with contours of the radio 
continuum at 4.86 GHz (total radio emission) observed with the 
VLA. The lowest contour is 50 /iJy beam"-*^ with a further spacing 
of 2 in ratio. Plot from Xu et al. (2003) reproduced by the friendly 
permission of the author and the A AS. 



outer tail could have formed earlier in time by an interac- 
tion event between NGC 7319 and an intruder, whereas the 
development of the inner tail seems to be caused by a more 
recent event (e.g. Moles et al. 1997; Xu et al. 2005), therefore 
the tails are sometimes also referred to the old and young 
tails. Note, however, that there are three different formation 
scenarios concerning the tails, which are still under discus- 
sion (see end of this section). The outer tail is extending 
to the southeast of NGC 7319 and is mostly covered by the 
foreground galaxy NGC 7320. The inner tail runs parallel to 
the outer tail below NGC 7319 and shows a higher surface 
brightness (Fedotov et al. 2011) and an active star formation 
(Xu et al. 2005). At its eastern end, the inner tail is pointing 
towards NGC 7320c. North of NGC 7318a and NGC 7318b 
there are two tidal arms extending from the galaxies, which 
belong to a region of active star formation (Gallagher et al. 
2001; Xu et al. 2005). A region of enhanced X-ray and radio 
emission is observed between NGC 7319 and NGC 7318b. 
This region is widely interpreted as a large shock front (cf. 
Fig. 1) resulting from a collision of the high-speed intruder 
NGC 7318b with the cold IGM (Shostak et al. 1984). 
Originally, the large shock was discovered in the 21 cm West- 
erborg radio continuum map by Allen & Hartsuiker (1972) 
(see Fig. 2). Later, the prominent ridge was shown to be also 
present on observed X-ray maps (Pietsch et al. 1997; Sulentic 
et al. 2001). Due to the additional, albeit faint, evidence of 
the ridge in observations within the optical wavelength, Xu 
et al. (2003) suggested that the shock front contains ionized 
gas and cold dust as well as hot thermal electrons, relativis- 
tic electrons and strong magnetic fields. However, the shock 
region between NGC 7319 and NGC 7318b is not the only 
radio source in SQ (cf. Fig. 2): further observed regions of 
enhanced radio emission are the innermost parts of NGC 
7319 and NGC 7318a. Also, the region north of NGC 7318a 
and NGC 7318b, which is observed to have a high star for- 



mation rate (see below), is also found to emit synchrotron 
radiation (Xu et al. 2003). 

The star formation activity in SQ is believed to be trig- 
gered by the interactions and can be measured in UV (Xu 
et al. 2005). The system reveals several regions of active star 
formation (^1 Mq yr~^), such as the inner disc of NGC 
7318b, the disc of NGC 7319 and the star-forming region 
in the north of the galaxy pair NGC 7318a/b. Cluver et al. 
(2010) found active star formation regions in the shock to 
be located in the intruder galaxy and at either end of the ra- 
dio shock, whereas the shock ridge itself seems to host only 
regions of weak star formation. Additional regions of weak 
star formation (^ 1 Mq yr~^) are found in the outer disc 
of NGC 7318b, within the inner and outer tails and in the 
overlap region of the tail (Xu et al. 2005). 
The outer and the inner tails are the result of interactions 
within SQ. Therefore, any model suggesting an interaction 
history for SQ should be able to reproduce the structure 
and age of these features. Three different possible interac- 
tion scenarios discussed in the literature meet this require- 
ment. The first scenario suggests a formation of the outer 
tail at an earlier time than the formation of the inner tail. 
Xu et al. (2005) proposed that the outer tail is the result 
of an encounter between the galaxies NGC 7320c and NGC 
7319. The inner tail is suggested to have formed later due 
to an encounter of NGC 7319 with NGC 7318a. The second 
possible scenario also suggests a formation of the tails at dif- 
ferent times but both resulting from passages of NGC 7320c 
over NGC 7319. Based on the fact that both the outer and 
the inner tails are directed towards NGC 7320c, Moles et 
al. (1998) and also Sulentic et al. (2001) proposed a bound 
orbit of NGC 7320c around NGC 7319, where the first inter- 
action created the outer and a subsequent passage the inner 
tail. The third scenario assumes a formation of both tails 
by only one interaction event of NGC 7320c and NGC 7319 
(e.g. Hwang et al. 2012) thus suggesting the same formation 
age for the two tails. Although the common opinion favors 
different formation times of the tails due to i.e. different star 
formation rates and cluster ages within the tails (e.g. Fedo- 
tov et al. 2011), this issue is still under discussion (Hwang 
et al. 2012). 



3 INITIAL CONDITIONS 

The first attempt towards a morphologically adequate rep- 
resentation of SQ in simulations was made by Renaud et 
al. (2010). They performed a large number of collisionless 
N-body simulations to find initial parameters for the four 
progenitor galaxies, including initial positions and velocities. 
Starting with these initial parameters, the system undergoes 
a number of interactions resulting in a morphological struc- 
ture comparable to observational findings. The large-scale 
configuration of the tidal features and the galaxies is gener- 
ally well represented. However, these simulations are purely 
gravitational and therefore not suitable for more detailed 
studies of e.g. intergalactic gas properties, shocks, star for- 
mation activity, magnetic fields, etc. Hence, Renaud et al. 
(2010) suggest to use their models as a basis for more com- 
plex simulations of SQ. 

Recently, Hwang et al. (2012) presented a further model, 
including also a gaseous component. They performed re- 
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Table 1. Parameters of initial galaxy setup for both models of SQ 



Galaxy Parameters 





NGG7319 


NGG7320C 


NGG7318a 


NGG7318b 


Concentration" CC 


12 


20/8 


12 


12 


Spin parameter A 


0.1 


0.15 


0.1 


0.1 


Disk mass fraction^ 


0.125 • Mtot 


0.05 • Mtot 


0.125 • Mtot 


0.125 • Mtot 


Bulge mass fraction^ mb 


0.0625 • Mtot 


0.0214 . Mtot 


0.0625 • Mtot 


0.0625 • Mtot 


Disk spin fraction 


0.125 


0.05 


0.125 


0.125 


Gas fraction f 


0.2 


0.2 


0.2 


0.2 


Disk height^ zq 


0.2 . la 


0.2 . la 


0.2 . la 


0.2 . la 


Bulge size^ 


0.2 . la 


0.2 . la 


0.2 . la 


0.2 . la 


Scale length of extended gas disc^ Ig 


6 • la 


6 • la 


6 • la 


6 • la 



SQ MODEL A 

Initial coordinates (x,y,z) [kpc/h] (0.0, 0.0, 0.0) (19.92, 10.45, 19.92) (-83.12, 0.0, 38.70) (-14.04, 17.64, -217.19) 

Initial velocities {vx,Vy,Vz) [km/s] (0.0, 0.0, 0.0) (-620.0, 232.5, -387.5) (465.0, -46.5, -93.0) (218.0, 0.0, 1025.0) 

Disc orientation (6, (/)) (0, 0) (0, 0) (180, 23) (0, -23) 



SQ MODEL B 

Initial coordinates (x,y,z) [kpc/h] (0.0, 0.0, 0.0) (8.9, -10.9, 10.9) (-50.0, 7.1, -40.0) (15.0, 1.4, -241.4) 

Initial velocities (1;^,%,^^) [km/s] (0.0, 0.0, 0.0) (35.9, 79.5, -77.5) (100.0, -27.0, -92.5) (20.0, -7.5, 350.0) 

Disc orientation (6, (/)) (0, 0) (0, 0) (180, 0) (180, 0) 



(a) Second column: model A / model B. 

(b) Mtot and l^j are given in Table 2. 



stricted three-body/SPH simulations of the SQ system using 
different models for the progenitor galaxies and a different 
formation scenario compared to Renaud et al. (2010). Yet, 
they were also able to reproduce the main tidal features of 
the system. Furthermore, they found indications supporting 
the hypothesis that the large shock between NGC 7319 and 
NGC 7318b has been caused by a high-speed collision of 
NGC 7318b and the IGM. They also studied the behaviour 
of gas clouds within the shocked region developing after 
the collision. They found a continuing production of small 
shocks in this region over a time span of ^ 10^ yrs. However, 
they were not able to study the star formation history and 
the gas temperature in more detail, because their treatment 
of heating and cooling was not accurate enough to draw 
reliable conclusions. Moreover, they did not perform any in- 
vestigations concerning X-ray emission, magnetic fields and 
radio emission. 

Both of the models include the four strongly interacting 
galaxies NGC 7319, NGC 7320c, NGC 7318a and NGC 
7318b. The galaxies NGC 7317 and NGC 7320 were not 
considered, as the influence of NGC 7317 on the other group 
members is negligible at present time (Renaud et al. 2010) 
and NGC 7320 is observed to be an unrelated foreground 
galaxy (Shostak et al. 1974; Allen & Hartsuiker 1980; Moles 
et al. 1997). 

Both Renaud et al. (2010) and Hwang et al. (2012) have 
considered the three different formation scenarios described 
in section 2. Renaud et al. (2010) found the flrst scenario to 
best represent the present-day morphology of SQ. They also 
found that the second scenario does not work for the sim- 
ulations at all. In contrast, Hwang et al. (2012) found best 



results using the third scenario. This is mainly because for 
the first scenario they were not able to reconstruct the tails 
one after another and to adjust the orbit of NGC 7318a in 
such a way that a good representation of the observations is 
reached. They ascribe the differences of the results mainly 
to the more extended halo potentials and also to the differ- 
ences between the purely gravitational N-body simulations 
and those including gas physics using SPH (Hwang et al. 
2012). 

Therefore, we chose in each case the best model, i.e. the first 
scenario model of Renaud et al. (2010) and the third scenario 
model of Hwang et al. (2012), and used them as a basis for 
SPMHD studies of SQ. The sequence of interactions within 
the model of Renaud et al. (2010) is the following: NGC 
7320c undergoes a collision with NGC 7319, then 7318a in- 
teracts with the already disturbed galaxy NGC 7319 and 
finally the high-speed intruder NGC 7318b hits the system. 
For the model of Hwang et al. (2012), the interaction his- 
tory is different: NGC 7320c performs a close passage around 
NGC 7319, then the galaxies NGC 7318a and NGC 7318b 
undergo a collision behind the orbital plane of the main sys- 
tem followed by a collision of the high-speed intruder NGC 
7318b with the IGM material west of NGC 7319. 

Modifications of the galaxy models and the initial positions 
and velocities were required due to the additional inclusion 
of a gaseous component and an ambient IGM. We note that 
we did not intend to reproduce the models of Renaud et 
al. (2010) and Hwang et al. (2012) in detail but rather to 
use them as basis for our studies of magnetic fields, X-ray 
and radio emission in SQ. In the following sections we pro- 
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Table 2. Galaxy Model Parameters. 



Model Mtot 
[IQIOMq] 


R200 
[kpc/h] 


[kpc/h] 


N 


-a 

halo 


^^disc 




^^bulge 






SQ 


MODEL 


A 


















NGC7319 


258.1 


199.0 


7, 


.1 


706 


900 


869 


600 


217 


400 


544 


000 


NGC7320C 


44.8 


111.0 


5, 


.9 


141 


380 


60 


952 


15 


238 


32 


610 


NGC7318a 


88.0 


139.0 


5, 


.0 


240 


346 


295 


660 


73 


920 


184 


960 


NGC7318b 


46.0 


112.0 


4, 


.0 


127 


242 


156 


530 


39 


130 


97 


920 






SQ 


MODEL 


B 


















NGC7319 


12.5 


63.0 


2, 


.6 


737 


456 


907 


636 


226 


908 


567 


272 


NGC7320C 


2.5 


42.5 


3, 


.6 


169 


780 


73 


136 


18 


284 


39 


128 


NGC7318a 


8.2 


63.0 


2, 


.3 


483 


364 


594 


908 


148 


728 


371 


820 


NGC7318b 


7.1 


60.0 


2. 


,2 


418 


008 


514 


472 


128 


620 


321 


544 



(a) collisionless particles within dark matter halo. (b) collisionless particles within disc, 
(c) gas particles within disc. (d) collisionless particles within bulge. 



vide details on our galaxy models and initial positions and 
velocities. 



3.1 Galaxy models 

Renaud et al. (2010) set up their galaxies to be composed 
of an exponential disc, a bulge and a dark matter halo (all 
consisting of collisionless gravitational N-body particles) us- 
ing a method based on Hernquist (1993). As NGC 7320c is 
assumed to be spherically symmetric it is only made up of 
a halo and a bulge component. 

For our first SQ model, we adopt the total masses as well as 
the percental masses of bulge, halo and disc (if present). The 
galaxies are set up using the method described by Springel 
et al. (2005b), which is also based on the Hernquist method 
(Hernquist 1993). This method allows for a galaxy model 
consisting of a cold dark matter halo, an exponential stellar 
disc, a stellar bulge (all of these components being colli- 
sionless gravitational N-body particles) and an exponential 
gaseous disc (SPH particles). We included a gaseous disc 
component using a disc gas fraction of f = 0.2 to all of the 
galaxy models. The number of particles and thus the reso- 
lution were highly increased. The parameters of our galaxy 
models are given in Table 1, whereby masses are given in 
units of the total galactic mass Mtot and scale lengths in 
units of the stellar disc scale length Id (Mtot and la are given 
in Table 2). 

The galaxy models of Hwang et al. (2012) are composed of 
a dark matter halo and a disc containing star as well as gas 
particles. NGC 7318a was indeed set up with a disc, but the 
angular and random velocities are representing an elliptical. 
Surprisingly, the total masses of the galaxies are roughly a 
factor of 10 smaller compared to the SQ model of Renaud 
et al. (2010). Also, the mass ratios of the galaxies among 
each other differ significantly. We discuss the effects of the 
smaller galaxy masses and different mass ratios within our 
simulations in chapter 6. 

For our second SQ model, we adopt the total masses of 
Hwang et al. (2012). However, most of the other parameters 



Table 3. Gravitational softening lengths e. 





Halo 


Disk 


Gas 


Bulge 




[pc/h] 


[pc/h] 


[pc/h] 


[pc/h] 


SQ model A 


113 


28 


28 


28 


SQ model B 


48 


10 


10 


10 



are chosen in accordance with Renaud et al. (2010). Again, 
our galaxy models consist of a cold dark matter halo, an 
exponential stellar disc, a stellar bulge and an exponential 
gaseous disc. Thus, compared to Hwang et al. (2012), our 
models additionally contain a bulge. We use a high resolu- 
tion of the galaxy models. As the true nature of NGC 7318a 
is still unclear, we have modeled it as a disc galaxy. Thus, 
the parameters of the galaxy models for our two SQ mod- 
els differ mainly in the mass (and size) of the galaxies. The 
parameters common to both of our SQ models are given in 
Table 1. 

For simplicity, we refer to the modified galaxy models with 
"galaxy model A" and accordingly with "SQ model A" to 
the total initial SQ setup including all galaxy models and an 
ambient IGM (see section 3.2) for our improvement based 
on the Renaud model and with "galaxy model B" and ac- 
cordingly "SQ model B" to our realization of the Hwang 
model. 

The total masses, virial radii and stellar disc scale lengths 
and particle numbers used in our galaxy models are sum- 
marized in Table 2. In case of model A, the setup results in 
particle masses of rrigas — ?Tidisc — mbuige ~ 2.1 • 10^ 
Mq and TTihaio ~ 2.1 • 10^ h-^ Mq with h = 0.7. For model 
B, the particle masses are mgas — '/^disc = '^^buige ~ 9-8 • 10^ 
Mq and mhaio ~ 9.8 • 10^ Mq. We note that 
these particle masses are of order of the largest molecu- 
lar clouds, i.e. small-scale shock turbulence is not mod- 
eled in our work. The fixed gravitational softening lengths 
(see e.g. Dehnen 2001; Johansson et al. 2009) are listed in 
Table 3. The softening lengths have been adjusted using 
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Table 4. Comparison of initial SQ models. 



SQ model of 
Renaud et al. (2010) 



SQ model A 



SQ model of 
Hwang et al. (2012) 



SQ model B 



DM halo 

Stellar bulge 

Stellar disc 

Gaseous disc 

Ambient gaseous IGM 

NGC 7318a modeled as elliptical 

Formation scenario 

Mass ratios of galaxies^ 

Order of total masses of galaxies 

Resolution: mass per star (gas)-particle 

Extension DM halo^'^ 



1 



: 0.2 : 0.3 : 0.2 

^ IO^^Mq 



! 3 • lO^Mf 
46 kpc 







/ 
/ 
/ 
/ 



: 0.2 : 0.3 : 0.2 

^ IO^^Mq 
=i 3 • IO^Mq 
280 kpc 



/ 

/ 
/ 

/ 
3 

: 0.2 : 0.7 : 0.6 

^ IO^OMq 

135 kpc 



/ 
/ 
/ 
/ 
/ 

3 

: 0.2 : 0.7 : 
^ 1010M( 
1 • IO^Mq 
102 kpc 



0.( 





(a) Mass ratio is given for NGC 7319 : NGC 7320c : NGC 7318a : NGC 7318b. 

(b) Cut-off radius for SQ model of Renaud et al. (2010) and Hwang et al. (2012), virial radius r200 for SQ model A and B. 

(c) The DM halo extensions are given for NGC 7319. 



enew = eoM-[{N^'VN''n • (MtX/Mt°oY)] The minimum 
SPH smoothing length for the gas particles is 1.0 e. 
Finally we want to note that our models differ from the 
correspondent models of Renaud or Hwang in some details: 
In our representation of the Renaud galaxy models (galaxy 
models A), the dark matter haloes of the galaxies are more 
extended and the disc scale lengths differ from the original 
model. These differences are on the one hand due to the 
enhancement of the original setup method for the galaxies 
(Hernquist 1993) by Springel et al. (2005b), and on the other 
hand due to the additionally included gaseous component. 
We feel confident that our new models are appropriate reali- 
sations of the SQ galaxies as they now contain gas and show 
more realistic extensions of the dark matter haloes. In our 
representation of the Hwang galaxy models (galaxy models 
B), the galactic disc sizes differ from the original model and 
we additionally include a bulge component. As in the re- 
stricted three-body simulations of Hwang et al. (2012) gas 
and star disc masses are negligible, we modeled the galaxies 
for our model B with parameters in analogy to galaxy mod- 
els A. The parameters which are common to both models 
are given in Table 1. A detailed comparison of the model 
properties is listed in Table 4. 



3.2 IGM 

In addition to the modifications described above, an ambient 
IGM is included. The IGM is set up to be composed of 
additional gas particles surrounding the galaxies similar to 
Kotarba et al. (2011). We arrange the IGM gas particles on a 
hexagonal closed-packed lattice. The particle masses of the 
IGM gas particles are adopted from the respective galaxy 
models. 

The IGM fills a volume of 1000 xlOOO xlOOO h'^ kpc^ and 
we assume a density of picM = 10~^^ g cm~^, resulting in 
particle numbers of A^igm= 496 828 for the SQ model A and 
iViGM= 2 311 367 for the SQ model B. 

The IGM is assumed to be already virialized, whereby the 
temperature within each model is set to the virial temper- 
ature at the virial radius of the largest galaxy model NGC 



7319: 

2 rripfi 1 2 .rripfx 

with the mean molecular weight for a fully ionized gas of pri- 
mordial composition /i ^ 0.588, proton mass nip and Boltz- 
mann constant ks- This leads to a temperature of the IGM 
of TiGM ~ 9.4-10^ K for the SQ model A and Tigm 1.2- 10^ 
K for the SQ model B. 



3.3 Initial positions and velocities of the galaxies 

Renaud et al. (2010) were the first to present possible N- 
body models of SQ. Hwang et al. (2012) resimulated their 
successful models including a gaseous component. However, 
their final model represents a different formation scenario 
than Renaud et al. (2010). 

Our models A and B differ in several features (e.g. gaseous 
component, IGM, dark matter distribution, see section 3.1) 
from the original models by Renaud et al. (2010) and Hwang 
et al. (2012), respectively. Therefore, we had to adjust the 
orbital parameters of the galaxies compared to the original 
models. We found these new orbital parameters of the galax- 
ies, i.e. new initial velocities in case of SQ model A and initial 
positions and velocities for SQ model B, by performing more 
than 100 test simulations to find out the best representation. 
The orbital parameters of the best representation of SQ for 
each model are found in Table 1. The parameters for the 
disc orientations are the same as in the original models. 



3.4 Initial magnetic fields 

The initial magnetization of the galactic discs is set up using 
= Bga\,o and By = Bz = G with the z-axis being the 
axis of rotation. This setup ensures that the initial field lies 
in the plane of the galactic disc. The initial magnetic field 
strength of the galaxies is assumed to be 5gai,o = 10~^ G. 
This value is by three orders of magnitude smaller than the 
typical observed galactic magnetic field value (e.g. Beck et 
al. 1996). For comparison, we also performed simulations 
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Table 5. Parameters for the multi phase model common to ah 
galaxy models. 

Multi Phase Model Parameters 



Gas consumption timescale 


tSFR 


8.4 Gyr 


Mass fraction of massive stars 




0.1 


Evaporation parameter 


Ao 


4000 


Effective supernova temperature 


TsN 


4 • 10^ K 


Temperature of cold clouds 


Tec 


1000 K 



with an initial galactic magnetic field strength of Bgai,o = 
IQ-^ G and Bgai,o = 10"^ G. 

The initial magnetic field of the IGM is assumed to be 
uniform in x-direction with an initial value of Bigm,o = 
i^iGM,x = 10~^ G and the x — y plane being the orbital 
plane. In this setup, the IGM magnetic field is naturally 
also pervading the magnetic field of the galaxies. As the 
intergalactic magnetic field is usually estimated to be less 
than 10~^ G, the assumed initial IGM magnetic field value 
is already close to the observed value. 



4 SIMULATION METHOD 

The simulations of SQ presented in this paper are performed 
with the N-body/SPH code Gadget (Springel et al. 2001; 
Springel 2005). Hydrodynamics is included with a formu- 
lation of SPH which conserves both energy and entropy 
(Springel & Hernquist 2002). The simulation results pre- 
sented in this work were obtained with the development 
version of Gadget-3, whereby the evolution of magnetic 
fields can be followed using the additional MHD implemen- 
tation of Dolag & Stasyszyn (2009). For a detailed descrip- 
tion of SPH and SPMHD methods see e.g. Dolag k Sta- 
syszyn (2009), Stasyszyn et al. (2012) and Price (2012). For 
all of the simulations, SPMHD is adopted using the standard 
(direct) magnetic field implementation, where the magnetic 
field is evolved using the induction equation. Additionally, 
within shocks the standard artificial viscosity and the artifi- 
cial magnetic dissipation are applied. For the dimensionless 
parameters we use values of av — 2.0 (viscosity), as — 0.5 
(dissipation) and /3 = 1.5 (e.g. Price 2012). We do not use 
a subgrid model for physical or turbulent magnetic dissipa- 
tion. 

The Lorentz force, describing the feedback on the plasma 
caused by the magnetic field, is taken into account as a con- 
tribution to the acceleration of each gas particle. However, 
in cases of strong magnetic forces this momentum conserv- 
ing form can lead to numerical instabilities, i.e. clumping 
of the particles (Phillips & Monaghan 1985), due to the 
non- vanishing numerical divergence of the magnetic field, 
which will not be stabilized when the magnetic pressure ex- 
ceeds the gas pressure (see e.g. Dolag & Stasyszyn 2009, 
Stasyszyn et al. 2012 and Price 2012). To avoid these numer- 
ical instabilities, the effects of any unphysical source terms 
of the magnetic field (i.e. any numerically non-vanishing di- 
vergence) are subtracted from the momentum equation as 
suggested by B0rve et al. (2001). A threshold for this diver- 
gence force subtraction (half the value of the current Lorentz 
force) is applied to account for situations in which the accel- 



SQ model A: 


Bgal = 10-9 




SQ model A: 


Bgal = 10-8 


G 


SQ model A: 


-Bgal = lO-'^ 




i 









SQ model B: Bgai = IQ-^ G 
SQ model B: Bgai = IQ-^ G 
SQ model B: Bgai = IQ-'^ G 




t [Myr] t [Myr] 

Figure 3. Mean numerical divergence {h\V • B|/B} as a function 
of time for (a) SQ model A and (b) SQ model B. Galactic mag- 
netic field (solid lines) and IGM magnetic field (dashed lines) are 
plotted separately using a threshold of 10"-^^ g cm""^. The nu- 
merical divergence measure stays essentially below the tolerance 
value of unity during all of the simulations. 



eration due to the divergence force could become dominant 
(see Stasyszyn et al. 2012). 

The MHD version of Gadget was already successfully em- 
ployed for studies of the magnetic field evolution in molec- 
ular clouds, in isolated and interacting galaxies as well as 
in galaxy formation simulations (e.g. Biirzle et al. 2011; Ko- 
tarba et al. 2009; Ceng et al. 2012; Beck et al. 2012). 
Radiative cooling, star formation (SF) and associated su- 
pernova feedback are applied using the hybrid multiphase 
model described in Springel & Hernquist (2003). The inter- 
stellar medium is modeled as a multi-phase gas consistent of 
condensed cold clouds embedded in an ambient hot gas at 
pressure equilibrium. These cold clouds are forming stars, 
of which a certain fraction is expected to die instantly as 
supernovae. Thereby, the feedback energy resulting from su- 
pernovae directly heats the ambient hot phase, which in turn 
looses energy by radiative cooling assuming a primordial gas 
composition and a zero-metallicity cooling function. The pa- 
rameters of the multi-phase model used for all of the galaxy 
models are chosen in accordance with Johansson et al. (2009) 
and are given in Table 5. For further details see the corre- 
sponding numerical papers on the Gadget code (Springel 
et al. 2001; Springel & Hernquist 2002; Springel 2005), the 
SPMHD method (e.g. Dolag k Stasyszyn 2009; Stasyszyn et 
al. 2012; Price 2012) or the star formation model (Springel 
k Hernquist 2003). 



5 SIMULATIONS 

Below, we present the results of our simulations of the mod- 
els A and B (cf. section 3). In all of our simulations, the 
numerical divergence {h\V • B|/B) (cf. Fig. 3) stays essen- 
tially below the tolerance value of unity (for a more detailed 
discussion of the numerical divergence, see e.g. Kotarba et 
al. 2010; Biirzle et al. 2011; Beck et al. 2012; Geng et al. 
2012), which implies the numerical reliability of the simula- 
tions. Unless specified otherwise, all the plots correspond to 
an initial magnetic field strength of Bigm = ^gai = 10~^ G. 
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Figure 4. (a) Gas density overlaid with contours of the total stellar surface density (contour levels: 0.003, 0.007, 0.03, 0.07, 0.3, 0.7, 3, 7 
M0 pc~^) for the present-day configuration of SQ model A (tA3 = 320 Myr). (b) Same as (a), but for the present-day configuration of 
SQ model B (f^s = 860 Myr) with three additional pink contour levels (levels: 0.00007, 0.0003, 0.0007 Mq pc~2), which were smoothed 
with a circular Gaussian beam with FWHM=3. (c) Gas density overlaid with the stellar surface density of new formed stars (contour 
levels: 3-10-'^, 0.003, 0.007, 0.02, 0.07 M© pc-^) for SQ model A (^^3 = 320 Myr). (d) Same as (c), but for SQ model B (^^3 = 860 
Myr). 



Table 6. Coordinates of the four galaxy models of SQ model A 
and SQ model B at the present-day configuration, respectively. 



SQ model A 


SQ model B 


(x,y,z) in [kpc/h] 


(x,y,z) in [kpc/h] 



NGC 7319 (-6.0, 3.0, 10.5) (5.0, 6.0, -17.0) 

NGC 7320c (-103.0, 14.0, -70.0) (-28.0, 30.0, -31.0) 

NGC 7318a (37.0, 11.5, -16.5) (30.0, -9.0, -84.5) 

NGC 7318b (29.0, 11.5, 35.0) (23.0, -3.5, -16.5) 



5.1 General morphology 

In the following, we denote particular evolutionary stages of 
our simulations with Al to A3 for the model A, and corre- 
spondingly Bl to B3 for model B in order to simplify later 
reference. 

For SQ model A, it takes 320 Myr from the initial config- 
uration (see Tables 1 and 2) to reach its best fit with ob- 



servations. In the course of the evolution, first NGC 7320c 
undergoes a collision with NGC 7319 (Al: tot = tAi = 80 
Myr), producing the outer tail. After tu = 140 Myr, NGC 
7318a starts interacting with the already disturbed galaxy 
NGC 7319, resulting in the formation of the inner tail. Af- 
terwards, the final encounter of the discs of NGC 7318a and 
NGC 7319 takes place at tA2 = 240 Myr. Subsequently, 
the high-speed intruder NGC 7318b hits the system, which 
leads about 40 Myr later to a configuration consistent with 
the morphology of the observed system (A3 - present day: 
tpd = tA3 = 320 Myr) (cf. Fig. 1). The resulting coordi- 
nates of the four galaxies at the present-day configuration 
are listed in Table 6. 

Fig. 4 (a) shows the gas density overlaid with contours 
of the stellar surface density for SQ model A at the time 
tA3 = 320 Myr (A3). Only NGC 7319 and the galaxy pair 
NGC 7318a/b are visible. NGC 7320c lies outside the plot- 
ted region (cf. Table 6). Qualitatively, the model shows a 
good agreement of the general features with observations, 
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Figure 5. SFR as a function of time for (a) SQ model A and (b) 
SQ model B. Note the different scaling of the axes. 



i.e. the morphology of the large galaxy NGC 7319, the in- 
ner tail south-east of NGC 7319, and the structure of the 
western galaxy pair NGC 7318a/b. Some kind of outflow is 
indicated in the gas density distribution north of NGC 7319 
and also north and south-west of the pair NGC 7318a/b. 
The stellar surface density shows the highest values within 
the inner disc of NGC 7319 and also within the disc of 
NGC 7318a. The stellar density in the disc of NGC 7318b is 
slightly lower. A bridge between NGC 7319 and the collid- 
ing pair NGC 7318a/b is clearly visible in the stellar surface 
density as well as in the gas density. At the left edge of 
the plot, there is an elongated region with a stellar surface 
density larger than 0.003 pc~^. This region belongs to 
the galaxy NGC 7320c which lies outside the plot. How- 
ever, the model does not reproduce the observed position of 
the galaxy pair NGC 7318a/b correctly and the outer tail 
is generated in this model but is already too diffuse to be 
visible at the present-day configuration. These aspects were 
both already noted for the original model by Renaud et al. 
(2010). Furthermore, the small-scale details of this galaxy 
pair, e.g. the structure of the spiral arms of NGC 7318b, are 
not well reproduced. 

For our SQ model B it takes 860 Myr from the initial condi- 
tions (cf. Tables 1 and 2) to develop a morphology similar to 
the observed configuration. However, the interaction history 
differs significantly from model A: At first, a close passage 
of NGC 7320c around NGC 7319 (Bl: tot = Ut = tsi = 160 
Myr) simultaneously forms both the inner and outer tails. 
Subsequently, the galaxies NGC 7318a and NGC 7318b un- 
dergo a collision behind the orbital plane of the main system 
consisting of the galaxies NGC 7319 and NGC 7320c (B2: 
tB2 = 640 Myr). Afterwards, the high-speed intruder NGC 
7318b is moving towards the main system and collides with 
the IGM material west of NGC 7319. Meanwhile, about 80 
h~^kpc behind the plane of the main system, NGC 7318a 
moves westwards. The evolution results in a configuration 
similar to observations (B3 - present day: tpd = tss = 860 
Myr). The resulting coordinates of the four galaxies at the 
present-day configuration are listed in Table 6. Within this 
model, NGC 7320c is still present in the tail, but however, 
it is largely disrupted, therefore we had to estimate the co- 
ordinates of the center of this galaxy. 

Fig. 4 (b) shows the same quantities as Fig. 4 (a) but for 
the SQ model B at the time tss = 860 Myr (B3). As for 
the model A, there is a good qualitative agreement of the 
general features with observations. Yet, model B reproduces 



the position of the galaxy pair NGC 7318a/b better than 
model A (cf. also Fig 1). The stellar surface density shows 
the highest values in the inner discs of NGC 7319 and NGC 
7318a/b. Again, a stellar bridge is visible between NGC 7319 
and the galaxy pair NGC 7318a/b. However, the position of 
NGC 7318a is slightly too much south, which is a feature 
of our SQ model B, not of the original model by Hwang 
et al. 2012. SQ model B is also not capable of reproducing 
small-scale features like the spiral arms of NGC 7318b or 
the detailed structure of NGC 7319 correctly and the outer 
tail is shorter compared to observations. 

SQ model A and SQ model B differ significantly in the for- 
mation scenario of the outer and the inner tails. Whereas 
in the first case the two tails are evolved in two different 
interactions of NGC 7319 with NGC 7320c and NGC 7318a, 
within the SQ model B the tails are created by only one 
interaction event of NGC 7319 with NGC 7320c. Within SQ 
model A, the outer tail is formed about tpd — tot = 240 Myr 
ago and the inner tail about tpd — tu = 180 Myr ago, re- 
sulting in a formation age difference of ^ 60 Myr (which is 
similar to the results of the model of Renaud et al. 2010, 
who found a formation age difference between the tails of 
^ 70 Myr). In contrast, SQ model B shows an equal for- 
mation age of both tails about tpd — tot = 700 Myr ago. 
However, observations are still suggesting different ages of 
the tails: Moles et al. (1997) proposed an age of the outer 
tail of ^ 500 - 700 Myr and for the inner tail - 200 Myr by 
considering the radial velocity difference between NGC 7319 
and NGC 7320c. Later, Sulentic et al. (2001) found that this 
measurement of the radial velocity was highly overestimated 
and suggested a much slower radial velocity for NGC 7320c 
(almost identical to that of NGC 7319) resulting in a predic- 
tion for the encounter of NGC 7319 with 7320c (causing the 
inner tail) within the second formation scenario (see chapter 
2) about ^ 500 Myr ago, which is similar to the predicted 
age of the outer tail. Fedotov et al. (2011) suggested an age 
of ^ 400 Myr for the outer tail and ^ 200 Myr for the inner 
tail. However, Fedotov et al. (2011) also found that the in- 
ner tail also contains some old clusters with an age of ^ 500 
Myr, even if it is mainly composed of blue clusters. Hwang 
et al. (2012) argue that the different formation ages would 
not have necessarily taken place in order to explain the dif- 
ferent star formation histories, as the outer tail evolves in 
a different environment, which may be less dense compared 
to the environment of the inner tail. 

SQ model B clearly supports the idea of a common origin 
of the both tails caused by only one interaction event about 
^ 700 Myr ago, which agrees well with observational pre- 
dictions for the formation age of the outer tail (e.g. Moles et 
al. 1997). Contrary, SQ model A implies different formation 
ages of the tails, whereby the finding of a formation age of 
the inner tail of ^ 200 Myr corresponds also to observational 
findings (e.g. Moles et al. 1997; Fedotov et al. 2011), whereas 
the formation age difference of only ^ 60 Myr seems to be 
too small compared to observations. To conclude, we cannot 
rule out any of the proposed formation scenarios and ages 
for the two tails, however, the slightly better agreement of 
the resulting positions of the galaxies (particularly of the 
galaxy pair NGC 7318a/b) and the larger formation age of 
the tails may be interpreted as an indication to prefer the 
formation scenario of SQ model B. 
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Figure 6. (a) Temperature overlaid with X-ray contours (logarithmic contour levels: 36.75, 37.25, 37.75 and 38.25 erg s"-*^) for the 
present-day configuration of SQ model A at the time t^s = 320 Myr. (b) Same as (a), but for the present-day configuration of SQ model 
B (ts3 = 860 Myr) and with lower contour levels (logarithmic contour levels: 35.75,36.25,36.75,37.25 erg s"-*^). (c) Magnetic field strength 
overlaid with arrows showing the direction of the magnetic field (initial magnetic field in x-direction with B^^\ = Biqm = 10~^ G) for 
SQ model A (t^s = 320 Myr). (d) Same as (c), but for SQ model B (Ibs = 860 Myr). 



5.2 Star formation 

Star formation (SF) activity in SQ is believed to be trig- 
gered by the interactions (Xu et al. 2005). Fig. 4 (c) again 
shows the gas density for model A at tAs = 320 Myr, but 
this time overlaid with contours of the stellar surface density 
of the new formed stars (i.e. stars which were not present in 
the initial setup). Most new formed stars are found within 
the inner disc regions of NGC 7319 and NGC 7318a. Slightly 
less SF takes place within NGC 7318b. A low star formation 
activity and thus a low surface density of newly formed stars 
is found within the spiral arms of NGC 7319, the outer discs 
of NGC 7318a/b, in the IGM between NGC 7319 and NGC 
7318a/b, and in the region of the inner tail south-east of 
NGC 7319. The surface density of new formed stars within 
the IGM mainly traces the stellar and gaseous bridge. How- 
ever, parts of this star- forming region seem to coincide also 
with probably the edges of a shock region visible in X-ray 
and radio emission (cf. section 5.3 and 5.5), which in prin- 
ciple corresponds to the observed star-forming regions (Xu 



et al. 2005; Cluver et al. 2010) (cf. section 2). However, the 
starburst region north of NGC 7318a/b found in observa- 
tions (Xu et al. 2005) is not revealed in the surface density 
of new formed stars in our SQ model A. 
Fig. 4 (d) shows the same quantities as Fig. 4 (c) but for 
the present-day configuration of SQ model B. The highest 
surface density of new formed stars is found within the inner 
discs of NGC 7319 and NGC 7318a/b. Less star formation 
takes place in the outer discs of the galaxies and within the 
spiral arm structure of NGC 7319. The lowest star forma- 
tion is found north of NGC 7319. Most of these regions are 
indeed also observed to form stars (Xu et al. 2005) (cf. sec- 
tion 2). However, there is no region of noticeable ongoing 
star formation in the IGM between the main galaxy NGC 
7319 and the galaxy pair NGC 7318a/b, i.e within or at the 
edges of the supposed shock region. The starburst region 
north of the pair NGC 7318a/b found in observations (Xu 
et al. 2005) is also not reproduced. 

The SF rate (SFR) of isolated galaxies is found to be approx- 
imately constant, whereby the constant SFR depends on the 
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Figure 7. Radial profiles of mean (solid lines), maximum (upper 
dotted lines) and minimum (lower dotted lines) temperatures for 
SQ model A at different evolutionary stages. Here, the origin is 
the centre of the main galaxy NGC 7319. 




Figure 8. Same as Fig. 7, but for SQ model B. The temperature 
gets enhanced by shocks and outflows, whereas for the more mas- 
sive galaxies of SQ model A this enhancement is more efficient. 



total masses and gas fractions of the progenitors (Cox 2004). 
This behaviour results in significantly higher starting values 
of the SFR in our SQ model A compared to SQ model B. 
Galaxy interactions trigger new SF as they efficiently com- 
press the gas and therefore lead to a conversion of cold gas 
into new stars (cf. section 4). The ability of the system to 
trigger prominent star bursts depends on the initial masses, 



mass ratios and on the orbit of the progenitor galaxies (Cox 
2004). Fig. 5 shows the total SFR as a function of the time 
for SQ model A (a) and SQ model B (b), respectively. In 
good agreement with previous studies of SFR during galaxy 
minor mergers (Cox 2004), the SFR strongly depends on 
the initial masses of the progenitor galaxies, resulting in a 
significantly higher SFR for SQ model A compared to SQ 
model B. In both simulations, the interactions significantly 
enhance the overall SFR, whereby the more massive galaxies 
of SQ model A lead to stronger star bursts compared to SQ 
model B. 

5.3 Temperature and X-ray emission 

The main source of X-ray emission is hot gas, which 
is heated by shocks accompanying the interactions. We 
calculate the bolometric X-ray luminosity following the 
method of Navarro et al. (1995), which assumes thermal 
bremsstrahlung to be the main X-ray source in agreement 
with the applied zero-metallicity cooling function. The bolo- 
metric X-ray luminosity is projected along the line of sight 
according to 
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(2) 

with mass mgas,i, density pi and temperature Ti of the i- 
th gas particle in cgs-units, respectively. Only fully ionized 
particles should be considered when calculating the lumi- 
nosity. Therefore, we exclude contributions of particles with 
temperatures lower than 10^'^ K and densities higher than 
O.OlMopc-^ (cf. Cox et al. 2006). 

Fig. 6 (a) shows the temperature overlaid with contours of 
the X-ray emission for the present-day configuration of SQ 
model A. The gas within the galaxies is cooler, whereas the 
IGM gas is heated by shocks and outflows caused by the in- 
teractions, which is in good agreement with previous studies 
(e.g. Kotarba et al. 2011; Ceng et al. 2012). The logarithmic 
contours illustrating the X-ray emission reveal a total X-ray 
luminosity which is approximately four orders of magnitude 
lower than the observed X-ray luminosity in SQ (10^° — 10^^ 
erg s~^, Sulentic et al. 2001). This low X-ray luminosity re- 
sults most probably from the applied zero-metallicity cooling 
and the lack of black holes in the simulations (cf. Cox et al. 
2006; Ceng et al. 2012). Nevertheless, the X-ray luminosity 
shows the highest values in the IGM region between NGC 
7319 and NGC 7318b, indicating a large shock east of NGC 
7318b. This shock region fits very well to the observed shock 
front visible as a ridge in the X-ray and radio emission (cf. 
section 2 and Fig. 2). However, the morphology of the shock 
region found in our simulations differs slightly from the ob- 
servations. This difference might be explained by the more 
northern position of the galaxy pair NGC 7318a/b compared 
to the observed position. 

Fig. 6 (b) shows the same quantities as Fig. 6 (a) for the 
present-day configuration of SQ model B. Again, the IGM 
surrounding the galaxies is heated by the shocks and out- 
flows caused by the interactions. One large outflow is clearly 
visible in the temperature in a region south of the galaxy pair 
NGC 7318a/b and a smaller outflow in the north-west of the 
pair. The gas within the galaxies is cooler than within the 
IGM. The overall temperature of this model is significantly 
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Figure 9. Mean total magnetic field strength as a function of time for (a) SQ model A and (b) SQ model B with initial magnetic fields 
of ^IGM = 10~^ G for each run and Bgai = 10~^ G (blue), Bgai = 10~^ G (green) and Bg^l = 10~^ G (red), respectively. Galactic 
magnetic field (solid lines) and IGM magnetic field (dashed lines) are plotted separately using a threshold of 10"-^^ g cm~^. In all of the 
simulations, the galactic magnetic field gets efficiently amplified, whereby the maximum magnetic field strengths during the interactions 
are higher for the more massive galaxies of SQ model A and for higher initial magnetic field strengths. 



lower compared to SQ model A (Figs. 6 (a) and 6 (b) use the 
same colour bar) . This is consistent with the smaller and less 
massive galaxy models. Consequently, also the X-ray emis- 
sion is approximately one order of magnitude smaller than 
for SQ model A. As before, the highest X-ray emission is 
found within the IGM region between NGC 7319 and NGC 
7318b, again indicating a prominent shock. This region of 
enhanced X-ray emission is more extended compared to SQ 
model A, and agrees in general well with observations (cf. 
Fig. 2). 

The radial profiles (with the origin being the centre of the 
main galaxy NGC 7319) of the mean, minimum and maxi- 
mum temperature within 260 Mpc at different evolutionary 
stages are shown in Fig. 7 for SQ model A and in Fig. 8 
for SQ model B. In both figures, the upper left panels show 
the profiles before any collision events at t = 20 Myr (AO 
and BO, respectively). Al, A2, A3 and Bl, B2, B3, respec- 
tively, denote the characteristic interaction events until the 
present-day configuration is reached (cf. section 5.1). For 
both models, at t =20 Myr (AO and BO), the temperature is 
lower at small radii, reflecting to the cooler gas within the 
galaxies. At greater radii, the mean temperature is constant. 
This is due to the initial setup, where the temperature was 
assumed to be already virialized (see section 3.2). Within 
the inner 50 - 100 kpc the small minimum values of 
the temperatures correspond to the gas within the smaller 
galaxies. In both models, the temperature gets enhanced by 
shocks and outflows caused by the interactions (A1-A3 and 
B1-B3, respectively). The increase in temperature is prop- 
agating towards higher radii due to the dilatation of the 
shock-heated regions (A2, A3 and B2, B3, respectively). The 
gas inside the galaxies is also heated by the interactions but 
it cools down again between the interactions. Consequently, 
the overall temperature inside the galaxies is generally much 
lower than the temperature of the IGM gas. As already de- 
scribed above, the larger and more massive galaxies of SQ 
model A enhance the IGM temperature more efficiently (Fig. 
7) compared to SQ model B (Fig. 8), in agreement with pre- 
vious studies (Geng et al. 2012). Moreover, the propagation 
of the shock-heated regions within the IGM is much more 



efficient for the SQ model A, because the smaller galaxies of 
SQ model B cause weaker shocks and thus a slower propa- 
gation of the shock- heated regions. 

5.4 Magnetic field structure 

The magnetic field is expected to get enhanced through ran- 
dom and turbulent motions driven by the interactions of the 
galaxies (see e.g. Brandenburg & Subramanian 2005, for a 
review). The compact SQ system has undergone a number 
of interactions, whereby the magnetic field should have been 
amplified significantly. Based on observational evidence, Xu 
et al. (2003) suggests that the shock front revealed by the X- 
ray and radio emission contains ionized gas and cold dust as 
well as hot thermal electrons, relativist ic electrons and mag- 
netic fields. This assumption agrees very well with expec- 
tations motivated by simulations: interaction-driven shocks 
are propagating favorably into the IGM, thereby heating the 
IGM gas and thus producing hot thermal electrons (Kotarba 
et al. 2011). Also, material of the galactic disc, i.e. ionized 
gas, is transported with the shock. Within the shock, elec- 
trons could become accelerated to relativistic velocities via 
the Fermi-acceleration process, in turn enhancing the radio 
emission in these regions. The IGM magnetic field is ampli- 
fied by shocks due to the compression of the field lines in 
front of the shock and turbulence behind the shock front. 
Fig. 6 (c) shows the mean total magnetic field strength over- 
laid with the magnetic field vectors for the SQ model A. It 
can be recognized that shocks and interaction-driven out- 
flows are expanding into the IGM, thereby enhancing the 
magnetic field strength up to values of almost /iC order. 
Some regions, such as the upper galactic arm of NGC 7319, 
the region around the inner tail, and parts within the galaxy 
pair NGC 7318a/b show lower values of the total magnetic 
field strength. These regions are most probably not directly 
affected by the encounters. Hence, the slight enhancement 
of the magnetic field in these regions might be solely due to 
the winding process of the galaxies. High values of the mag- 
netic field strength are in particular found within the large 
outflow in the upper north which seems to originate between 
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Figure 10. Radial profiles of mean (solid lines), maximum (upper 
dotted lines) and minimum (lower dotted lines) magnetic field 
strengths (with the origin being the centre of the main galaxy 
NGC 7319) for SQ model A at different evolutionary stages. 
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Figure 11. Same as Fig. 10, but for SQ model B. In both models, 
the magnetic field strength within the galaxies gets efiiciently am- 
plified during the interactions. The amplification is more efficient 
for SQ model A. 



NGC 7319 and NGC 7318b, where also the large shock is 
located. There, we also find high magnetic field strengths. 
The direction of the magnetic field, which is heading towards 
NGC 7318b, may indicate that the shock was triggered by 
a collision between NGC 7318b with the IGM (instead of 
a collision with NGC 7318a), thus supporting the common 
hypothesis of the origin of the shock (cf. section 2). 
Fig. 6 (d) shows the same quantities as Fig. 6 (c) but for the 



SQ model B. Note the different color bar appropriate for 
the lower magnetic field strengths. As before, shocks and 
interaction-driven outflows are propagating into the IGM, 
enhancing the magnetic field strength up to values of or- 
der 0.1 /iC. Again, regions of lower magnetic field strength 
are found between the tidal arms of NGC 7319 and within 
the two galaxies NGC 7318a/b. These regions are proba- 
bly not directly affected by the interactions (the galaxies 
NGC 7318a/b are affected by their mutual collision only 
in the very outer regions of the discs). The highest values 
of the magnetic field strength are found between the main 
galaxy and the galaxy pair NGC 7318a/b. Again, this highly 
magnetized region coincides with the supposed shock region 
also observed in X-ray. Comparable magnetic field strengths 
are found within the disc of NGC 7319 and within the out- 
flows north of NGC 7319 and south-west of NGC 7318a/b. 
The magnetic field vectors within the region of the shock 
are either directed towards NGC 7318b or along the inner 
tail. As the collision of NGC 7318b with NGC 7318a hap- 
pened about 220 Myr before the present-day configuration, 
the prominent shock found in SQ model B is certainly a 
result of an interaction of NGC 7318b with the IGM. 
Fig. 9 shows the evolution of the mean total magnetic field 
strength as a function of time for the SQ model A (a) and 
B (b) with initial magnetic fields of Bigm = 10~^ G within 
each run and Bga\ = 10~^ G (blue lines), Bg^i = 10~^ G 
(green lines) and Bgai = 10~^ G (red lines), respectively. We 
plot the galactic magnetic field (solid lines) and the IGM 
magnetic field (dashed lines) separately using a threshold 
of 10~^^ g cm~^. The presented simulations show a slight 
amplification of the total galactic magnetic field before the 
first encounter due to the winding process, which is in good 
agreement with previous studies (e.g. Kotarba et al. 2009, 
2010). During the first interaction of NGC 7319 with NGC 
7320c {tAi = 80 Myr and Ibi = 160 Myr, respectively) the 
galactic magnetic field gets efficiently amplified, whereby the 
maximum magnetic field strength reached during this inter- 
action are higher for the more massive progenitor galaxies 
of SQ model A. Also, they are higher for higher initial mag- 
netic field strengths, in agreement with previous simulations 
(Ceng et al. 2012). Within all of the simulations, further in- 
teractions of the galaxies (A2 and A3, B2 and B3, cf. section 
5.1) lead to further increase of the galactic magnetic field. 
This further amplification is again much more efficient for 
the more massive galaxy models of SQ model A. At the time 
of the present-day configuration, the galactic magnetic field 
strengths still show an increasing trend and thus have not 
yet reached the saturation levels (dynamic equilibrium or 
equipartition between turbulent and magnetic energy den- 
sity) of the galactic magnetic field within the different mod- 
els. We note that the total galactic magnetic field is a mean 
value of all four galaxies which had undergone interactions 
of different intensity and violence. The IGM magnetic field 
of both SQ models gets only slightly amplified during the 
encounters and reaches an value of order 10~^ to 10~^G, 
depending on the simulation. Again, the interactions of the 
more massive SQ galaxy models A lead to a more efficient 
amplification of the IGM magnetic field. 
The radial profiles (with the origin being the centre of the 
main galaxy NGC 7319) of the mean, minimum and max- 
imum magnetic field strength within a radius of 260 Mpc 
are shown in Fig. 10 for SQ model A and in Fig. 11 for SQ 
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Figure 12. Radio maps showing the total Intensity hot of model A (left panels) and model B (right panels), respectively, at 4.86 GHz 
(upper panels), 1.4 GHz (middle panels) and 0.2 GHz (lower panels). Colors visualize the gas density (color coding same as in Fig. 4a 
(left) and Fig. 4b (right). Upper panels: (a) the contour levels are 50, 100, 200, 400, 800, 1600 /xJy/beam and (b) 0.5, 1, 2, 4, 8, 16 
/iJy/beam. Middle panels: (c) the contour levels are 90, 180, 360, 720, 1440, 2880 /iJy/beam and (d) 1.2, 2.5, 5, 10, 20, 40 /iJy/beam. 
Lower panels: (e) the contour levels are 300, 600, 1200, 2400, 4800, 9600 /iJy/beam and (f) 2.5, 5, 10, 20, 40, 80 /xJy/beam. The pink 
contour in each plot corresponds to 50 /iJy/beam. We assume a beamsize of 20'^ Within both models, the pink contour level encloses 
larger areas for lower frequencies. 
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Figure 13. Synthetic radio map for model A (a) and model B (b), respectively, at 4.86 GHz. Colors visualize the total intensity Itot (in 
^iJy/beam). Black contours show the polarized Intensity Ip and correspond to (a) to 25, 50, 100, 200, 400, 800, 1600 /iJy/beam and (b) 
to 0.25, 0.5, 1, 2, 4, 8, 16 /iJy/beam. Magnetic field lines derived from calculations of polarization are shown in white. Note the different 
color coding of the radio maps. Within both models, high total and polarized synchrotron emission is found within the region of the 
prominent shock between NGC 7319 and NGC 7318a/b. 



model B, respectively, at different evolutionary stages. Both 
simulations show an amplification of the magnetic field at 
small radii, i.e. in the region where the galaxies are located, 
prior to the first interaction (AO and BO, respectively). This 
amplification is due to the winding process. At larger radii, 
the magnetic field strengths stay constant at the value of 
the initial IGM magnetic field (10~^G). In both models, 
the magnetic field strength at smaller radii, i.e. within the 
galaxies, gets efficiently amplified during the interactions. 
As before, for the more massive SQ model A, the amplifica- 
tion is more efficient. Also, within SQ model A, the higher 
magnetic field is transported farther out into the IGM due 
to the stronger outflows and thus a stronger dilatation of 
the shocked regions. 

5.5 Radio emission and polarization maps 

We compute synthetic radio emission and polarization maps 
for our simulations of the two different representations of SQ 
in an analogous way as presented in Kotarba et al. (2011). 
The calculation is outlined in Appendix A. 
Fig. 12 shows the synthetic radio maps where contours of the 
total intensity are overlaid on a density color plot (we use 
the same color code as in Fig. 4a and Fig. 4b, respectively). 
The left panels (a, c, e) show the present-day configuration 
of the SQ model A, and the right panels (b, d, f ) the present- 
day configuration of the SQ model B. The radio emission is 
calculated for three different frequencies = 4.86 • 10^ Hz, 
zy = 1.4 . 10^ Hz and zy = 0.2 • 10^ Hz (from top to bottom). 
For the SQ model A, the contour levels of the radio emission 
are similar to the levels used by Xu et al. (2003) (cf. also 
Fig. 2). As SQ model B shows a much lower radio emission, 
the contour levels are different. In all panels, we highlight 
the 50 /iJy/beam level with a pink contour. 
The synthetic radio maps of SQ model A (left panels in Fig. 
12) reveal a high total synchrotron intensity in the central 
part of NGC 7319, in the disc of NGC 7318a and in the re- 



gion between NGC 7319 and the galaxy pair NGC 7318a/b. 
The region of enhanced synchrotron intensity east of NGC 
7318b coincides very well with the shock region found in the 
observed X-ray luminosity. The regions of high synchrotron 
emission correspond well to observations (cf. Fig. 2, and Xu 
et al. 2003, their Fig. 4). There is a further region of en- 
hanced synchrotron emission north of NGC 7318b, which 
is presumably due to a strong outflow. The strength of the 
synchrotron intensity corresponds to the observed order of 
magnitude, but however, a direct comparison of the abso- 
lute strength of the synchrotron intensity is rather difficult 
due to the simplifying assumption of the cosmic ray (CR) 
energy distribution (cf. Appendix A). 

The synthetic radio maps of SQ model B (right panels of 
Fig. 12) reveal a high total synchrotron intensity within the 
disc of NGC 7319 and a slightly lower synchrotron emis- 
sion around NGC 7318a. Also, there is a region of high syn- 
chrotron emission between NGC 7319 and NGC 7318a/b, 
which coincides with the shock region found in the X-ray 
emission. In general, the regions of enhanced synchrotron 
emission agree well with observations (cf. Fig. 2, and Xu 
et al. 2003, their Fig. 4). However, the strength of the syn- 
chrotron emission within this SQ model is approximately a 
factor of 100 lower compared to SQ model A. There is a re- 
gion of slightly enhanced synchrotron intensity in the north 
of the galaxy pair NGC 7318a/b, which may be compared 
to the observed star-forming region also showing high radio 
emission (Xu et al. 2003). 

Within both models, the contour level of 50 /iJy/beam 
(highlighted in pink) encloses larger areas for lower fre- 
quencies. This behaviour is caused by two cumulative ef- 
fects: First, the assumed power-law spectrum of the CR 
electrons (see Appendix A) results in fewer CR electrons 
with high energies compared to lower energies. The CR 
electrons with high energies are responsible for the high- 
frequency synchrotron emission, whereas the CR electrons 
with lower energies lead to a low-frequency emission. Conse- 
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Figure 14. Radial profiles (with the origin being the centre of the 
main galaxy NGC 7319) of mean (solid lines), maximum (upper 
dotted lines) and minimum (lower dotted lines) synchrotron emis- 
sion at 1/ = 4.86 • 10^ Hz for SQ model A at different evolutionary 
stages. 




Figure 15. Same as Fig. 14, but for SQ model B. Within both 
models, the total synchrotron intensity gets amplified during the 
interactions. The amplification is more efficient for the more mas- 
sive SQ model A. 



quently, the low frequency emission should be more intense 
(more photons) than the high frequency emission. Second, 
the strongest magnetic fields develop in the region of the 
large shock, the galaxies and the outflows. Thus, the high- 
frequency synchrotron emission should be concentrated in 
these areas. Together, this leads to an overall low intensity 



of the high-frequency emission. On the other hand, there are 
wide areas of weak magnetic field producing a high amount 
of low- frequency emission. Thus, the low-frequency emission 
is at last more intense. Applying smaller frequencies, we find 
a slightly higher visibility of details in the synchrotron dis- 
tribution. But generally the regions of enhanced synchrotron 
emission remain comparable to the observations of Xu et al. 
(2003, cf. their Fig. 4) (cf. Fig. 2). 

Fig. 13 shows the synthetic radio map of model A (a) and 
model B (b), respectively, at 4.86 GHz for the present-day 
configuration (A3 and B3, respectively, cf. section 5.1). Col- 
ors visualize the total intensity Itot and the black contours 
show the polarized Intensity Ip (in /iJy/beam). The direc- 
tion of the magnetic field is indicated by the white lines, 
whereby the length of these lines is scaled according to the 
degree of polarization Hobs (cf. length scale in lower left 
corner of the plots). Total as well as polarized synchrotron 
emission is visible inside but also outside the galactic discs, 
indicating that the interaction-driven shocks have already 
magnetized the sourrounding IGM. In both models, high 
total and polarized synchrotron emission is found within 
the region of the large shock between NGC 7319 and NGC 
7318a/b. The contour levels of the polarized synchrotron 
emission within SQ model B are approximately a factor of 
100 lower compared to SQ model A. 

The radial profiles (with the origin being the centre of the 
main galaxy NGC 7319) of the mean, minimum and maxi- 
mum total synchrotron emission at zy = 4.86 • 10^ Hz within 
a radius of 260 Mpc are shown in Fig. 14 for SQ model 
A and in Fig. 15 for SQ model B, respectively, at different 
evolutionary stages. The profiles of the synchrotron inten- 
sity evolve similar to the profiles of the magnetic field (Figs. 
10 and 11): Within both models there is an amplification 
of the total synchrotron emission at small radii prior to the 
first interaction (AO and BO). At larger radii, the total syn- 
chrotron emission stays constant, which corresponds to the 
constant initial IGM magnetic field strength and orientation 
in the outer regions. During the interactions, the total syn- 
chrotron intensity gets amplified within both SQ models. At 
the time of the second and third interaction (A2 and A3, B2 
and B3), the synchrotron emission is already entering larger 
radii, i.e. the IGM. In case of the more massive SQ model A, 
the initial synchrotron intensity within the galaxies and the 
IGM is higher and also the amplification is more eflficient 
compared to SQ model B. Also, the synchrotron emission is 
transported farther outwards into the IGM. 



6 DISCUSSION OF THE SQ MODELS 

The comparison of our numerical simulations with observa- 
tions of SQ yields some agreement with several observed 
properties of the gaseous and stellar components of SQ. 
However, the quality of the match with observations is dif- 
ferent for our two SQ models. In Table 7 we present a brief 
listing of the achievement of our simulations in comparison 
with the results of the previous numerical studies of SQ by 
Renaud et al. (2010) and Hwang et al. (2012). 
In general, both presented models of SQ are capable of repro- 
ducing its large-scale structure, but the spacial distribution 
of matter is different due to the different sizes of the partic- 
ipating galaxies. Furthermore, the enhancement of gaseous 
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Table 7. Achievement of our SQ simulations compared to previous numerical studies of SQ 





SQ model by 
Renaud et al. (2010) 


SQ model A 


SQ model by 
Hwang et al. (2012) 


SQ model B 


SQ large-scale morphology 


qualitative agreement 
with observations 


qualitative agreement 
with observations 


qualitative agreement 
with observations 


qualitative agreement 
with observations 


Presence of large shock 




clearly visible in X-ray 
and radio emission 


indicated by gas distribution 
and shock-heated gas 


clearly visible in X-ray 
and radio emission 


Star forming regions 




general agreement 
with observations 


hardly comparable 
(very roughly displayed) 


correspond roughly 
with observations 


Star formation rates 




behaviour in agreement with 
previous numerical studies 




behaviour in agreement with 
previous numerical studies 


Gas temperature 




behaviour in agreement with 
previous numerical studies 


hardly comparable 
(very roughly displayed) 


behaviour in agreement with 
previous numerical studies 


X-ray emission regions 




correspond generally 
with observations 




correspond generally 
well with observations 


X-ray emission strengths 




too low compared to 
observations 




too low compared to 
observations 


Magnetic field 




behaviour in agreement with 
previous numerical studies 




behaviour in agreement with 
previous numerical studies 


Radio emission regions 




good agreement with 
observations 




correspond generally to 
observations 


Radio emission strengths 




order of magnitude corresponds 
to observations 




low compared to 
observations 



properties is naturally more efficient for SQ model A com- 
pared to SQ model B, which is mainly due to the larger total 
masses of the involved progenitor galaxies. 

Regions of high synchrotron intensity might correspond to 
regions of high density, however, within periods of intensive 
shock ejection, the magnetic field is amplified by turbulence 
and shocks. During the interactions, the magnetic field and 
correspondingly the total and polarized synchrotron inten- 
sities are enhanced behind the shocks and transported into 
the IGM. This behaviour is well displayed in our SQ simu- 
lations. 

The observed prominent shock within SQ develops within 
both of our SQ models. The shock region is found in both the 
synthetic X-ray and the radio maps of our simulations. The 
distribution of the total synchrotron intensity compares very 
favorably with observations (cf. Figs. 2 and 12, see also Xu et 
al. 2003). The highest synchrotron values are reached within 
the region of the large shock in both SQ models. However, 
model A reproduces the strength of the synchrotron emission 
in the observed order of magnitude, whereas the extension 
of the shock region is slightly better reproduced by model B, 
which clearly underestimates the strength of the synchrotron 
emission. 

The final configuration in SQ model A is reached after the 
high-speed intruder NGC 7318b hits the system and finally 
collides with NGC 7318a. Nevertheless, we claim that the 
large shock within this system results from the collision of 
NGC 7318b with the IGM. Especially the direction of the 
magnetic field vectors (which are directing towards NGC 
7318b) supports this origin of the shock, which is also con- 



sistent with general acceptance. However, it cannot be ruled 
out that the large shock in SQ model A results at least par- 
tially from a collision of NGC 7318b with NGC 7318a. On 
the contrary, the large shock in SQ model B is certainly a 
result of an interaction of NGC 7318b with the IGM as the 
interaction between NGC 7318a and NGC 7318b happened 
approximately 220 Myr before the present-day configuration 
is reached. 



Overall, we can recognize a general trend in our simulations: 
the enhancement and propagation of regions with higher 
values of the studied gaseous properties (i.e. temperature. 
X-ray emission, magnetic field strength and synchrotron in- 
tensity) resulting from interactions and associated shocks 
and outflows is much more efficient for the more massive 
galaxies of SQ model A. SQ model B shows qualitatively 
similar effects of enhancement and propagation, but to a 
smaller extent. This is not surprising, as higher masses lead 
to higher kinetic energies of the interacting galaxies, which 
results in higher equipartition values i.e. the thermal en- 
ergy (temperature) or the magnetic energy (magnetic field 
strength) . These general findings and the good agreement of 
the synchrotron intensity in model A with observations may 
be interpreted as an indication for larger progenitor galaxies 
within SQ. Therefore, the promising extension of the shock 
region in both X-ray and synchrotron emission motivating 
the underlying formation scenario of SQ model B in com- 
bination with the total masses of the SQ model A would 
provide a very good starting basis for further studies. 
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7 CONCLUSIONS 

We have presented simulations of Stephan's Quintet includ- 
ing magnetic fields, radiative cooling, star formation and su- 
pernova feedback. We have investigated different properties 
of the gaseous component for two different galaxy models 
based on Renaud et al. (2010), SQ model A, and based on 
Hwang et al. (2012), SQ model B, respectively. We have set 
the focus on the general morphology, on the distribution of 
star-forming regions and star formation rates, on the tem- 
perature and the corresponding X-ray emission and finally 
on magnetic fields and the resulting total and polarized radio 
emission. A brief listing of the achievement of our simula- 
tions in comparison with the previous studies by Renaud 
et al. (2010) and Hwang et al. (2012) is shown in Table 7. 
The main results of our simulations can be summarized as 
follows: 

• The present-day configuration of SQ model A devel- 
ops within 320 Myr. The morphology of the system agrees 
qualitatively well with observations, only the position of the 
galaxy pair NGC 7318a/b, its small-scale details and the 
inner and outer tails cannot be reproduced correctly. The 
outer tail is generated in this model but already too dif- 
fuse to be visible at the present-day configuration as already 
noted for the original model of Renaud et al. (2010). 

• The present-day configuration of SQ model B develops 
within 860 Myr. Again, the morphology of the system agrees 
qualitatively well with observations, however, the position 
of NGC 7318a is slightly too southern. Also, the small-scale 
features such as the arms of NGC 7318b or the smaller-scale 
structure of NGC 7319 cannot be reproduced correctly and 
the outer tail is shorter compared to observations. 

• Within SQ model A, the total masses of the galaxies 
are approximately Milky Way- like. In contrast, the galactic 
masses of SQ model B are roughly a factor of 10 smaller 
compared to SQ model A. As lower galactic masses imply 
lower equipartition energies, the enhancement of the gaseous 
properties is commonly lower for SQ model B. 

• The regions of active star formation within SQ model A 
are found mainly in the discs of the galaxies, and also within 
the inner tail and between NGC 7319 and the pair NGC 
7318a/b. The latter partly coincides with the region of the 
large shock. Within SQ model B, the regions of active star 
formation are found within the inner discs of NGC 7319 and 
the galaxies NGC 7318a/b, but there is no region of active 
star formation between these galaxies. 

• The global SFR strongly depends on the initial masses 
of the galaxies and is significantly higher for the more mas- 
sive SQ model A. 

• In both models, the temperature of the gas within the 
galaxies is cooler compared to the IGM, which gets heated 
by shocks and outflows caused by the interactions. The mean 
temperature in SQ model B is significantly lower compared 
to SQ model A. 

• The X-ray emission shows the highest luminosities in 
the region of the large shock between NGC 7319 and the 
pair NGC 7318a/b within both models, in good agreement 
with observations (Pietsch et al. 1997; Sulentic et al. 2001). 
The X-ray luminosity in the shock region within SQ model 
B is about one order of magnitude smaller compared to SQ 
model A. 

• We find high values of the magnetic field strength in the 



region of the large shock and also within outflow regions in 
both SQ models. The values of the magnetic fleld strength 
within SQ model B are approximately a factor of 3 smaller 
compared to SQ model A. 

• The temporal evolution of the mean total magnetic fleld 
reveals an ampliflcation of the magnetic fleld strengths along 
with the different interactions. Thereby, the increase in the 
magnetic fleld strength is more efficient for SQ model A. 

• The synthetic radio maps of both models show a high 
total and polarized synchrotron intensity within the large 
shock, within NGC 7319 and around and within NGC 7318a. 
This finding agrees well with observations (cf. Xu et al. 
2003). 

The large shock revealed by observations of SQ is most likely 
the result of a colhsion of NGC 7318b with the IGM. The 
observed ridge of radio emission can therefore be ascribed to 
shock activity. The shock front in our simulations is clearly 
visible in the X-ray and synchrotron emission within both 
SQ models. We emphasize the importance of shocks for the 
magnetic field amplification and the enhancement of the syn- 
chrotron emission. Whenever a high amount of synchrotron 
emission is detected in regions between interacting galaxies, 
it may be ascribed to shock activity. 

For future studies, a further development of the existing 
SQ models would be essential to draw more detailed con- 
clusions on the extension and strength of the synchrotron 
emission within SQ. As the SQ model B results in a lower 
enhancement of the gaseous properties mainly because of the 
smaller masses, but displays the regions of enhanced X-ray 
and synchrotron emission quite well, it would be worthwile 
to use a different scaling of the total masses of SQ model 
B comparable to the total masses of the SQ model A. This 
would lead to a better comparability of the strengths of the 
gaseous properties of the present-day configuration of the 
two different models of SQ. Another particular focus in fur- 
ther studies should thereby be placed on the position and 
extension of the galaxy pair NGC 7318a/b, which we found 
to significantly affect the extension and structure of the large 
shock in SQ. As in our simulations the used particle masses 
are of the order of the mass of the largest molecular clouds, 
small-scale turbulence within the large shock region as re- 
cently observed by Guillard et al. (2012) cannot be modeled 
in our work. Therefore, further numerical simulations focus- 
ing on smaller scales would lead to a deeper understanding 
of the involved processes of shock activity, especially shocks 
wrapped around clouds and cloud like structures. Further- 
more, observations of the radio emission ridge at different 
frequencies would be of particular interest in order to gain 
new insights into the shock region. This knowledge could 
then be used as a basis for further improvements of numer- 
ical SQ models. 
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APPENDIX A: CALCULATION OF TOTAL 
AND POLARIZED SYNCHROTRON EMISSION 

The calculations of the total and polarized synchrotron emis- 
sion were performed using the code P-Smac2 (Donnert et 
al., in preparation), which projects the corresponding val- 
ues on a grid. For convenience of the reader, the required 
calculations are outlined below. 

The radiation spectrum can be evaluated for a given dis- 
tribution of electron energies. The energy spectrum of the 
cosmic ray electrons can be approximated by a power-law: 

n{E)dE = KE-^'dE, (Al) 
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with n{E) being the number density of electrons and p = 2.6 
the index of the power spectrum. The latter implies a spec- 
tral index a = {p — l)/2 of 0.8, which was also found by 
(Xu et al. 2003) for NGC 7319. is a constant normaliza- 
tion factor of the cosmic ray energy spectrum, which can be 
calculated for a given CR energy density ecR- 



ecR 



-f 



En{e)dE 



E^^-^'^dE, (A2) 



with an assumed energy range of Emin = 10^ eV to Emax = 
10^^ eV. The CR energy distribution is expected to be pro- 
portional to the thermal energy distribution. We assume 

eCR = 0.01 • etherm- 

We calculate the total synchrotron emission Jjy at a given 
frequency v via (Longair 2011) 



xT 



TleC^ip + 1) 
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with the Gamma function F, the electron mass me, the elec- 
tron charge e and the speed of light c. The frequency is set 
to u = 4.86 • 10^ Hz and u = Hz, respectively, which 

correspond to the frequencies used for the observations pre- 
sented by Xu et al. (2003) (see Fig. 2) . Additionally, we also 
use a frequency of i/ = 0.2 • 10^ Hz. B± denotes the mag- 
netic field component perpendicular to the line of sight and 
is taken from the simulations. 

Within each grid cell, the total synchrotron emission Ji, is 
calculated. Integration of Ji, along the line-of-sight gives the 
total synchrotron intensity Itot- The polarized emission Ipoi 
can be obtained using the Stokes parameters Q and U : 



^pol 



Jloi 



Jiy cos(2'0)ds 



Jloi 



J^sin(2^)ds ,(A4) 



where the integration is carried out over the line-of-sight 
(los). In case of a power- law form of the energy spectrum, 
the degree of polarization can be evaluated using (Longair 
2011) 



H 



p±l_ 



(A5) 



The polarization angle ip is defined as the angle between 
the electric field vector E ± of the radiation perpendicular 
to the magnetic field and the x-axis in the plane of the sky 
(xy-plane). Thus, sm(2ip) and cos(2'0) can be calculated as 
follows: 



sin(2^) : 



2Bx By 



COS(2'0) : 



Bl + Bl' 



(A6) 



The observed degree of polarization Hobs can finally be cal- 
culated using: 

-pr _ Ipol 

J-tot 

To avoid errors due to the spatial isotropy of the emission, 
the total and polarized intensities are multiplied with a fac- 
tor 



(A7) 



with d being the distance to the observer and rbeam the as- 
sumed radius of the beam corresponding to a resolution of 
20" . For the distance to SQ, we use the estimated value of 
94 Mpc (Moles et al. 1998; Appleton et al. 2006). The artifi- 
cial flux then corresponds to the expected amount reaching 
Earth from the distance of SQ. 
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